Hermitian Schmidt Decomposition and Twin 
Observables of Bipartite Mixed States 



F Herbutf 

Faculty of Physics, University of Belgrade, POB 368, Belgrade 11001, Yugoslavia and 
The Serbian Academy of Sciences and Arts, Knez Mihajlova 35, 11000 Belgrade 

Abstract. Study of mixed-state quantum correlations in terms of opposite- 
subsystem observables the measurement of one of which amounts to the same as that of 
the other, of so-called twins, is continued. Twin events that imply biorthogonal mixing 
of states, called "strong twin events", are studied. It is shown that for each mixed 
state there exists a Schmidt (super state vector) decomposition in terms of Hermitian 
operators, and that it can be the continuation of the mentioned biorthogonal mixing 
due to strong twins. The case of weak twins and nonhermitian Schmidt decomposition 
is also investigated. For separable states a necessary and sufficient condition for 
the existence of nontrivial twins is derived. Utilization of the Hermitian Schmidt 
decomposition for finding all twins is illustrated in full detail in the case of the two 
spin-one-half-particle states with maximally disordered subsystems (mixtures of Bell 
states). It is shown that only rank two mixtures have nontrivial twins. 



PACS numbers: 3.65.Bz, 03.67.-a, 03.67.Hk 



Submitted to: Institute of Physics Publishing 



f E-mail: fedorh@infosky.net 



Hermitian Schmidt Decomposition and Twins of Mixed States 



2 



1. Introduction 

Nowadays one distinguishes sharply between separable bipartite mixtures, which are 
quasiclassically correlated, and nonseparable ones, endowed with entanglement, a purely 
quantum property. (A good example of the latter is the case of correlated pure 
states.) The term "quantum correlations" is used in the generic sense comprising both 
quasiclassical correlations and entanglement. 

It was claimed in a recent investigation [1] that the study of quantum correlations 
through twin observables, or shortly twins, is expected to be important for quantum 
communication and quantum information theories because it is believed to reveal 
some basic properties of the correlations. Twin observables are opposite-subsystem 
observables such that the (subsystem) measurement of one of them amounts ipso facto 
to a measurement also of the other. Equivalently put, the subsystem measurement of a 
twin gives rise, on account of the quantum correlations, to an orthogonal decomposition 
of the state of the opposite subsystem. 

In bipartite mixed states it is easier to relate twins to quantum correlations than 
to entanglement (though the latter is more important). A quantitative measure of the 
former is, what is called, von Neumann's mutual information 

C(pi 2 ) = S(p 12 \pi ® p 2 ) = S(pi) + S(p 2 ) - S(p 12 ), 

expressed in terms of the so-called relative (or conditional) entropy, and, alternatively, 
in terms of the von Neumann entropies of the reduced statistical operators pi, i = 1, 2 
(subsystem states) and the von Neumann entropy of the statistical operator (bipartite 
state) pi2 itself. I denote it by "C" because it was thus designated and called 
"logarithmic correlation" by Lindblad 2j. He also made use of the classical discrete 
mutual information I(A, B\pi 2 ) of two arbitrary opposite subsystem observables A and 
B (with purely discrete spectra) which were assumed to be simultaneously measured in 
a quantum state p\ 2 . Then, utilizing the spectral forms and the ensuing probabilities 

A = Y,^P[ k \ k^k' => a h ±a v , B = J2biQ { 2 \ => h^hy- 

k I 

p(k,l) = Tr[p 12 (Pl k) ®Q 2 l) )}, p k = J2p(Kl), Pi = ^2p(k,l); 

l k 

one defines the mutual information 

H(A : B\p l2 ) = H(p(k,l)\p kPl ) = H( Pk ) + H( Pl )-H(p(kJ)), 

where H(p k ) e. g. is the so-called Gibbs-Boltzmann-Shannon entropy H(p k ) = 
— J2 k Pklogp k , etc. Finally, Lindblad defined 

I(A,B\p l2 ) = supH(A:B\p l2 ), 

where the supremum was taken over all possible choices of the observables. 
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Lindblad showed that 

I(A,B\p l2 )<C(p 12 ), 

and that 

C( Pl2 )>0 I(A,B\p 12 )>0 

are always valid. 

Thus, in all correlated states, i. e., in states in which C(p\ 2 ) > 0, or equivalently, 
P12 7^ Pi ® P2j one can understand part of the quantum correlations in terms of 
simultaneous subsystem measurements and their maximal mutual information. 

Now, twins take up a very special position among the subsystem observables, 
because, if A and B are twins then 

I(A,B\p 12 ) = H(p k ), 

since H{p k ) = H{pi) = H(p(k, I)) due to p(k, I) = Pk^ij{k)i where f{k) is a fixed bijection 
of the values of k onto those of I. This is the case of perfect correlations, called "lossless 
and noiseless information channel" in information theory. 

The investigation of twins began with pure states [3] jlj p\ 2 =|$)($|. Surprisingly, a 
necessary and sufficient condition for a subsystem observable A to have a nontrivial twin 
was found in terms of properties of p\ alone (local properties, cf (2a)). The opposite 
subsystem observable B that is the twin of A was, naturally, expressed in terms of global 
properties of | $). These were in a simple way given in terms of an operator (called 
correlation operator cf (11)) mapping the range of p\ onto that of p 2 . It was defined 
by | $). This operator is most practically handled in terms of the so-called Schmidt 
decomposition [3] because it is precisely the (antiunitary) operator determining which 
characteristic vector of p 2 should appear in the same term as that of p\ in the mentioned 
decomposition [3] (cf (11)). 

When twins were investigated in the mixed-state case [I], the mentioned condition 
(cf (2a)) was found to be only necessary. Actually, a sufficient condition for A to have 
a twin expressed as a property of p\ alone (a local property) cannot exist because for 
every p\ there is the uncorrelated state p\ 2 = p\ <8> P2, which does not have nontrivial 
twins. 

Thus, global properties inherent in pi 2 have to be made use of at the mentioned very 
first stage of investigation of twins in the mixed state case. It is not easy to "extract" 
a minimal global property of p\ 2 that does the job (as in the pure state case). 

It is a striking fact that the Schmidt decomposition of state vectors can be 
generalized to all mixed states. It is the basic aim of this article to investigate the 
relevance of this decomposition to twins. It is proved that the Schmidt decomposition 
of any bipartite mixed state p\ 2 need not be expressed in terms of some very general 
linear operators, it can be given exclusively in terms of Hermitian operators, which can, 
in principle, be physically interpreted as observables (cf Theorem 2 and Corollary 1). 

The concept of strong twins, which are closely connected with biorthogonal 
decomposition of p\ 2 (cf Theorem 1), is introduced as a step towards the mentioned 
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Hermitian Schmidt decomposition of p\ 2 . Also nonhermitian Schmidt decomposition of 
mixed states is studied (cf Theorem 3). 

For mixtures of Bell states a Hermitian Schmidt decomposition is given in the 
literature (though not treated as such). On this simple example the problem of finding 
all twins is easily solved (cf Theorems 5 and 6) in order to illustrate the relevance of the 
Hermitian Schmidt decomposition to extracting the sought for global property inherent 
in Pi2- 

In the mentioned simple case it turns out that rank four mixtures do not allow 
nontrivial twins. This is not surprising because it was shown in the preceding study 
[T] that singularity of pn is a necessary condition. But, surprisingly, also rank three 
mixtures are shown to have no nontrivial twins. This suggests that perhaps a stronger 
necessary condition, some kind of "sufficient singularity", for the existence of nontrivial 
twins could be found in the general case. This will be followed up elsewhere. 

Relating twins to separability is fully clarified in this study in terms of a necessary 
and sufficient condition for the existence of nontrivial twins (cf Theorem 4). Relating 
twins to entanglement in the mixed-state case, and to the quantitative measures of 
entanglement like the so-called entanglement of creation and entanglement of distillation 
jHj, or the quantum relative entropy [7] and others is an important open question that 
will be hopefully treated in further work. 

The study of twins having been pursued in a number of mentioned articles is an 
ab ovo approach, which is already proved to be, in principle, relevant and perhaps 
even important to quantum information theory. It stands somewhat apart from the 
mainstream investigations. But it will be, hopefully, connected up with the latter as a 
result of further exploration. 

2. Preliminary Relations 

When a general, i. e., mixed or pure, bipartite state (statistical operator) p i2 is given, 
twins (Ai, A 2 ) are algebraically defined as Hermitian (opposite subsystem) operators 
satisfying 

Mf>X2 = A 2 p 12 , (1) 

where A\ is actually (Ai <S> h), h being the identity operator for the second subsystem, 
etc. It was shown pQ that (1) implies 

[A 1 ,p 1 } = 0, [A 2 ,p 2 ] = (2a, b) 

for the subsystem states (the reduced statistical operators). (The symbols Tr i? i = 1, 2, 
denote the partial traces.) Relation (2a) is the mentioned local necessary condition on 
Ai to have a twin. 

If Pi is a first-subsystem projector, one can decompose the statistical operator: 



Pl2 = P\P\2 + A Pl2, 



(3) 



Hermitian Schmidt Decomposition and Twins of Mixed States 



5 



where Pj 1 is the orthocomplementary projector of P\. Let (Pi, P2) be a pair of nontrivial 
twin events (twin projectors) for p 12 . In general, the terms on the RHS are not even 
Hermitian. First, we are going to investigate the more important case when (3) is a 
mixture of states. 

3. Strong twin projectors and biorthogonal mixtures 

Let (Pi,P 2 ) be a pair of nontrivial twin projectors for a composite-system statistical 
operator p 12 - 

Remark 1. Evidently, either both terms on the RHS of (3) are Hermitian or none of them 
is. They are Hermitian if and only if the projector P 1 (or equivalently, Pj 1 ) commutes 
with p i2 : 



(any one of the equalities implies the other), as seen by adjoining the terms in (3). 

Hermiticity of the terms in (3) implies that they are statistical operators (up to 
normalization constants), i. e., that (3) is a mixture. Namely, if (4) is valid, then 
idempotency leads to P\pn = P1P12P1, which is evidently a positive operator. Since 



the operator has a finite trace. 

Definition 1. Nontrivial twin events (projectors) we call either strong twin events (pro- 
jectors), if they satisfy (4), or weak twin events (projectors) if (4) is not satisfied. 

A strong twin event P\ implies a mixture (3) of states that have a strong property 
called biorthogonality. To understand it, we first remind of (ordinary) orthogonality of 
states. 

If p' and p" are statistical operators with Q' and Q" as their respective range 
projectors, then one has the known equivalences: 



[P,Pi 2 ] = 0, 



1,2, 



(4) 



TrP^Px < Trp 12 = 1, 



p'p" = Q'Q" = <£> TZ(p') JL K{p"), 



(5) 



where the last relation expresses orthogonality of the ranges. 

Any of the three relations in (5) defines orthogonality of states. 



Definition 2. If 



P12 = wp' 12 + (1 - w)p\ 



'l2) 



< w < 1, 



(6) 



is a mixture of states such that 



p' iP '! = 0, i = l,2, 



(7) 
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where p[ = Tr 2 p' 12 etc. are the reduced statistical operators, then we say that (6) is a 
biorthogonal mixture. 

To prove a close connection between strong twin events and biorthogonal mixtures, 
we need another known general property of composite-system statistical operators p i2 : 

P12 = Q1P12 = P12Q1 = Q2P12 = P12Q2, (8) 

where Qi is the range projector of the corresponding reduced statistical operator p iy 
i = 1,2. 

Theorem 1. If P x is a nontrivial twin event, (3) is a biorthogonal mixture if and only if 
Pi is a strong twin event. 

Proof. Sufficiency. If P\ is a strong twin projector and (6) is obtained by rewriting (3), 
then p' 12 = Pip' 12 is valid, and this implies p[ = P\p\ for the reduced statistical operator, 
and, adjoining this, one arrives at p[ = p' x P\. On the other hand, one has analogously 
P12 = P \P\2 implying p'[ = P^p'[. Finally, 

p'ip'i , = (p , iA)(P 1 Vi , ) = o. 

The symmetrical argument holds for the second tensor factor. 

Necessity. If (6) is a biorthogonal mixture, then we define Pj = Q[, i = 1,2, i. e., 
we take the range projectors of the reduced statistical operators of p' 12 as candidates for 
our twin projectors. On account of (8), we can write (6) as follows: 

P12 = wQWMQ'i + (1 - w)Q'[Qy[ 2 Q'[Qi. 

Since in view of (5) biorthogonality (7) implies Q'iQ" = 0, i = 1, 2, it is now obvious 
that Pi and P 2 , multiplying from the left pi 2 , give one and the same operator, i. e., that 
they are twins, and it is also obvious that they both give the same irrespectively if they 
multiply pi 2 from the left or from the right, i. e., that they are strong twin projectors. 

□ 

In view of (5), it is clear that biorthogonal decomposition of a statistical operator 
can be, in principle, continued: If, e. g., p' 12 in a biorthogonal decomposition (6) is, in 
its turn, decomposed into biorthogonal statistical operators and replaced in (6), then 
any two of the new terms are biorthogonal etc. 

An extreme case of a biorthogonal mixture is a separable one: 

Pl2 = E^(Pl fc) ®P2 fc) ), (9) 
k 

where 

VA; : w k > 0, pf } > 0, Trpf } = 1, 1 = 1,2; ^> fc = 1 

k 
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("p > 0" denotes positivity of the operator). This decomposition cannot, of course, 
always be carried out, but examples are well known. For instance, if one performs ideal 
nonselective measurement of the z-component of spin of the first particle in a singlet 
two-particle state, one ends up with 

p 12 = (l/2)(|z+) 1 (z+| 1 ® \z-) 2 (z-\ 2 + \z-) 1 (z-\ 1 ®\z+) 2 (z+\ 2 ). 

This is obviously a biorthogonal separable mixture. 

One wonders if, at the price of relaxing the requirement of statistical-operator terms 
as slightly as possible, there could exist a general decomposition into uncorrelated terms 
(like in (9)). 

To find an affirmative answer, we take resort to the known case of general (entangled 
or disentangled) composite-system state vectors and their Schmidt decompositions. Let 
us sum up the sufficiently detailed relevant information on this [3]. 

The Schmidt decomposition of an arbitrary pure state vector | $)i 2 of a composite 
system is expressed in terms of its canonical entities. They are the following: 

(i) The reduced statistical operators (subsystem states) pi ( = Tr 2 |$)i2($|i2 ) and 
p 2 (defined symmetrically) are well known. 

(ii) The spectral forms of the reduced statistical operators are 

Pi = J2 r i l*)i(*li> P2 = J2 r i N>2<* U, Vi : n > 0. (10a, 6) 

i i 

(Note that the positive spectra -multiplicities included - are always equal.) 

(iii) Finally, the mentioned expansion utilizes the (antiunitary ) correlation operator 
U a , which maps the range 7Z(pi) onto the range lZ(p 2 ). (Note that they are always 
equally dimensional in the pure state case). The correlation operator is determined by 
|$)i2, and, in turn, in conjunction with pi, it determines |$)i2- 

The Schmidt decomposition reads: 

l<f)i2 = E^ 1/2 N)i®(^N)i) 2 - (ii) 

i 

The normalized characteristic vectors \i) 2 in (10b) may (and need not) be chosen to be 
equal to (ll a I *)i) 2 - 

In case of a state vector |$)i 2 , the characteristic relation (1) for twins reduces to 

Ai |$)i2 = A 2 |$) 12 . (12) 

The corresponding twin A 2 then satisfies 

A 2 = U a A 1 U- 1 Q 2 + A 2 Q^, (13) 

where Q 2 is the range projector of p 2 , and Q2, its orthocomplementary projector, 
projects onto the null space of p 2 . 

One should note that, on account of the commutation (2b), both the range and 
the null space of p 2 are invariant for A 2 . Further, the second term on the RHS of (13), 
or rather the restriction of A 2 to the null space, which corresponds to it, is completely 
arbitrary and immaterial for the twin property (12), because it acts as zero on | $)i 2 . 
(Naturally, the symmetric claim holds true for A\ and pi) 
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4. Hermitian Schmidt decomposition of bipartite statistical operators 

It is well known that linear Hilbert-Schmidt operators A acting in a Hilbert space, i. 
e., those with a finite Hilbert-Schmidt norm (TrAU) , form a Hilbert space in their 
turn. Writing the operator A as a (Hilbert-Schmidt) supervector \A), the scalar product 
is 

(A\\B}= TiA^B. 

Since for every statistical operator p, one has Trp 2 < 1, it is a Hilbert-Schmidt 
operator. Therefore, every statistical operator has a Schmidt decomposition. 

The trouble is that the operators that take the place of the state-vector tensor 
factors in the terms of (11), which are the sought for generalizations of the statistical 
operators p\ , i — 1, 2 in (9), are in general linear operators. This might be a too wide 
generalization. One wonders if one could be confined to Hermitian operators. 

When we view the operators as supervectors, then we must view adjoining of 
operators as an antiunitary operator the square of which is the identity operator, i. e., 
which is an involution. Hence, we denote adjoining by 1 ®V 2 a ^ for a composite system. 
Hermitian are the operators that are invariant under the action of this antiunitary 
involution. 

Fortunately, the Schmidt decomposition can always be expressed in terms of Her- 
mitian operators. We put this in a more precise and a more detailed way. But it is 
simpler to return to the Hilbert space of state vectors to perform some elaboration. 

Theorem 2. Let v} a) ®vj a) be a given antiunitary involution acting on composite-system 
state vectors. One has the equivalence: 

(V} a) ® V 2 {a) ) | $} 12 =| $> 12 & [ Pu Vf a) ] = 0, • = 1, 2; V 2 {a) U a Vl a) = U a , (14) 

where pi, U a are the above mentioned canonical entities of | $)i 2 . (Note that in the last 
relation we, actually, have the restriction of Vi to IZ(pi).) 

Proof. Let |$)i 2 be invariant under the action of the antiunitary involution. Then 

V^pM^ = n (a) (Tr 2 |$>i 2 ($h )vl a) = 

Tr 2 (vl a) |$) 12 ($| 12 V{ a) ) = Tr 2 {v} a) [(v} a) ® vf } ) |$) 12 ($| 12 (v} a) ® v} a) )] v} a) } = 
Tr 2 (v 2 {a) |$)i 2 ($|i 2 y 2 (a) ) =Tr 2 | $) 12 ($ | 12 = p u 



and symmetrically for p 2 . One has to note that an antiunitary involution equals its 
inverse and its adjoint. Further, use has been made of some known basic properties of 
partial traces (which are analogous to the well known ones for ordinary traces). 

Commutation of p\ with the antiunitary involution allows one to choose the 
characteristic basis {| i)i : Vi} of the former spanning its range consisting of vectors 
invariant under the action of (cf jHj). 
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Now, let us take the Schmidt decomposition (11) in terms of such an invariant 
basis. Then 

(Vt ] ® Vi a) ) | $)i2 = £r? /2 l<>i ® V 2 {a) (u a |z)i) 



Since |$)i2 is assumed to be invariant, it follows that also 

l*>ia = E''* 1/2 |» , >i®V 2 (B) (cr 8 |i> 1 ) a . 

The second tensor factor in each term is uniquely determined by the LHS and the 
corresponding first tensor factor (as a partial scalar product, cf [3 ). Comparison with 
(11) then shows that 

Vi: V 2 {a) U a \i) l = U a \i) l . 
Since = Vi \i)i, we further have 

V (") UaV f«) = Ua 

as claimed. 

Conversely, if the main canonical entities are in the relation to the antiunitary in- 
volutions as stated in (14), then we can expand | $)i 2 in a characteristic basis of p\ 
spanning its range that is invariant under the antilinear operator. Then (11) immedi- 
ately reveals that, as a consequence, | $}i 2 is invariant under Vi ® V 2 . □ 

Corollary 1. Every composite-system statistical operator p\ 2 has, after normalization, 
a Hermitian Schmidt decomposition. 

Proof. Since every pi 2 , being Hermitian, is invariant under the antiunitary involution 
(g> V 2 \ Theorem 2 immediately implies that p\ 2 , upon super vector normalization, 
has a Schmidt decomposition in terms of Hermitian operators. □ 

Returning to a biorthogonal mixture, one wonders if one can continue such a 
decomposition by writing each term in a Hermitian Schmidt decomposition in order 
to obtain the latter decomposition for the entire statistical operator. The answer is 
affirmative on account of the following: 

Going back to (5), we can add a fourth equivalent property. 

Proposition 1. Two statistical operators p' and p" are orthogonal if and only if they 
are orthogonal as Hilbert- Schmidt supervectors. 



Proof. It is obvious that orthogonality (in the sense of (5)) implies Hilbert-Schmidt 
orthogonality. To see the converse implication, we make use of the fact that every 
statistical operator has a purely discrete spectrum |HJ, and we decompose the statistical 
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operators in terms of characteristic vectors corresponding to positive characteristic 
values: 

(p' 1 1 P") = Tl p'p" = Tl J2 r k\ k )( k \J2 f i \j)(j\ = 

k j 

k j 

Hence, 

(p'||p">=0 p ' p " = 

(cf the third relation in (5)). □ 

If (A\, A 2 ) is a pair of twin observables, then the detectable parts A' { , i = 1,2, have 
a common purely discrete spectrum {a n : Vn} (with, in general, different multiplicities), 
and the corresponding (detectable) characteristic projectors {P} : % = 1,2 Vn}, are 
also pairs of twins pQ. 

Definition 3. If a// mentioned characteristic projector pairs (P[ , P^) are strong twin 
projectors, then (Ai,A 2 ) is a pair of strong twin observables. If some of the detectable 
characteristic twin projectors are strong and some weak, we say that we have partially 
strong (or, synonymously, partially weak) twin observables. If all the mentioned twin 
projectors are weak, then we have a weak pair of twin observables. 

A pair (Ax, A 2 ) of nontrivial twin observables for p 12 is a pair of strong ones if and 
only if 

[A i ,p 13 ] = 0, * = 1,2 (15) 

is valid. This is so because commutation with all characteristic projectors is equivalent 
to commutation with the Hermitian operator itself, and, if Pi e. g. is a nondetectable 
characteristic projector of A\, then one has commutation because 

P1P12 = (PiQi)pi2 = = PniQiPi) = PuPi 

on account of (8). 

Strong twin observables, by means of their strong characteristic twin projectors, 
lead to a generalization of (3): 

Pl2 = E P i n Vl2 = J2 W nPl2, (16a) 
n n 

where 

\fn: w n ^Tip 12 P[ n \ p^ = (w n r 1 P^p 12 . (166) 

Naturally, if P[ P12 = 0, then p ™ 2 is not defined. Any two terms in (16a) are 
biorthogonal. 

Note that we utilize the entire characteristic projectors, which are the orthogonal 
sums of the detectable and the nondetectable parts: P[ n ^ = (P[)^ © (P")^ parallelling 
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Ux = K( Pl ) © K L (pi) because (P x ') (n) Pi2 = Pi (n W 
Proposition 2. If 

Jn) _ T (n) 

Pi = tr 2 pi 2 5 

and symmetrically for p 2 , are the reduced statistical operators of the terms in a 
biorthogonal mixture (16a), then 

Pt ] pt ] =pt\ < = 1,2, (17a) 

or equivalently, 

7^(pS n) )C7e(^ (n) ), 1 = 1,2. (176) 

Proof. On account of the definition of (16a), one has Pj^ p%2 = P12 ■ Taking the 
opposite-subsystem partial trace, one obtains P^pf 1 ' = p^ i = 1,2. □ 

Corollary 2. If the detectable part A[ of a twin observable A\ has a nondegenerate 
characteristic value a n corresponding to a strong characteristic twin projector 
(PQO) — I ^p( n )^ 1 (^ip( n ) | 1; | t//™)^ g IZ(pi), then the term in the biorthogonal mixture 
(16a) that corresponds to it has the form 

WniV^W^Ilrf, (18) 

where p 2 n ^ is a (second-subsystem) statistical operator and (18) is a term in a final Her- 
mitian Schmidt decomposition of pi 2 . 

Any biorthogonal decomposition of a composite-system statistical operator p 12 
(into two or more terms) can be continued in each term separately into a Schmidt 
decomposition of P12 in terms of Hermitian operators. 

The biorthogonal decomposition is an intermediate step. This is similar to the case 
when we can partially diagonalize the Hamiltonian of a quantum system (due to some 
symmetry e. g.). The diagonalization is then continued separately with each submatrix 
on the diagonal of the Hamiltonian. 

The continuation from a biorthogonal mixture to a Hermitian Schmidt 
decomposition can always be performed, in principle, "by brute force": diagonalizing 
the reduced statistical superoperator Pl of the normalized supervector |p 12 ) (analogously 
as it is done for an ordinary state vector), and by finding an invariant basis for in 
each characteristic subspace thus obtained |Hj. 

5. Weak twins and nonhermitian Schmidt decomposition 

For the sake of completeness it is desirable to investigate decomposition (3) also for 
a weak nontrivial twin projector P\. First, we take an analytical view of Theorem 1 
to realize that the biorthogonality of the two terms in (3) is connected with the twin 
property (strong or weak), and the strong twin property corresponds to the hermiticity 
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of the terms. Let us put this more precisely. 
Remark 2. A decomposition 

P12 = A 12 + B X2 

of a composite-system statistical operator p 12 into two linear operators is biorthogonal 
if there exist two opposite-subsystem projectors (Pi,P 2 ) such that 

A 12 = PiA 12 = P 2 A 12 , = P 1 B 12 = P 2 B 12 ; 

= PfA 12 = P£A 12 , B 12 = PfB 12 = P 2 L B 12 . 

It is clear from Theorem 1 that any birthogonal mixture (of states) (6) satisfies the 
condition given in Remark 2. Having in mind (3), it is also evident that biorthogonality 
is equivalent to the existence of a pair of twin projectors (weak or strong). Finally, the 
strongness property of the twins is equivalent to the hermiticity of the terms in (3), 
which results in having statistical operator terms (and a mixture). 

Theorem 3. If (Pi, P 2 ) is a pair of weak twin projectors for a composite-system statistical 
operator pi 2 , then the terms in (3) are super vectors, and replacing each by a (nonher- 
mitian) Schmidt decomposition, one obtains a decomposition of the same kind for the 
entire statistical operator. 



Proof. Since in 



1 > Trp^ 2 = Trpi 2 Pipi2 + Tip 12 P 1 p 12 



the terms are nonnegative (as traces of positive operators), the terms in (3) are Hilbert- 
Schmidt operators, i. e., super vectors. Suppose we have decomposed the first term in 
(3) in the Schmidt way: 

P lP n = cj:rl /2 A? 



B 2 , 



where c is a normalization constant (because the statistical operator is not a super state 
vector unless it is a pure state). Since the LHS is invariant under P 1 , so is each first- 
subsystem linear operator A± \ because the second factors in the expansion have unique 
corresponding first factors. If we decompose also the second term in (3) in the Schmidt 

way 



AV2 = C'EW 



u 2 i 



then, analogously, invariance of each factor under P^~ follows. This results in super 
vector orthogonality: 



Vi, j : 



Tr 



(A?yc[ j) \ =Tr (A { i ) yP 1 Pj L C { 1 



= 0. 



The symmetrical argument goes for the second factors and P 2 . Thus, replacing both 
terms in (3) by their nonhermitian Schmidt decompositions, we have biorthogonality 
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between any term of the first decomposition and any term of the second one. Therefore, 
we have a decomposition of the same kind of the entire pw a 

It is now clear that also in the case of weak twin projectors the decomposition (3) 
can be continued, but this time to a nonhermitian Schmidt decomposition. 

A nonhermitian Schmidt decomposition need not be wild and far fetched from 
the physical point of view. Let me illustrate this by the obvious fact that a Schmidt 
decomposition of a state vector | $)i 2 

l$>i2 = E r ' /2 l*>2> V *^' : (*l P |i% = 0, P=l,2 

i 

immediately results in a nonhermitian Schmidt decomposition of the statistical operator 
|$>i2($|i 2 : 

l$>i 2 ($li2=EE r i /2r v /2 l*)i<*'li ® l*>2<*' U • 

i i' 

Finally, let us return to separable mixtures. 
6. Nontrivial twin projectors for separable mixtures 

Let (9) be a general separable mixture. Let us clarify under what conditions it has 
nontrivial twin events. 

Theorem 4- A general separable mixture (9) has a nontrivial twin projector P\ if and 
only if the set of all values of the index "k v is the union of two nonoverlapping subsets, 
say, consisting of "A/" values and of "k"" values respectively, and, when (9) is rewritten 
accordingly: 

Pl2 = E W k'p { l } ® P2 } + E W k"Pl } ® P2 \ (190) 
k' k" 

then one has biorthogonality between the two groups of terms: 

W, W: pPpP = 0, z = 1,2. (196) 
Before we prove the theorem, we first prove subsidiary results. 

Lemma 1. Let 

Pl2 = E W -l Vl/(m) )l2( V I /(m) |l2 

rn 

be an arbitrary pure-state mixture. Then, a pair of subsystem observables (Ai, A 2 ) are 
twins for p±2 if and only if they are twins for all pure-state terms. 

Proof. Necessity follows from the general result that all twins of pu are also twins of all 
state vectors from the topological closure TZ^pu) of the range of pyi (cf section 3, CI in 
lj). As well known, the vectors {|^/ ( ' Tn ^)i2 : Vm} span the mentioned subspace. 
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Sufficiency is obvious. □ 
Lemma 2. Let 

k 

be an arbitrary mixture. The pair (Ai,A 2 ) are twin observables for p\ 2 if and only if 
they are twin observables for all term states pi 2 ■ 

Proof is immediately obtained from Lemma 1 if one rewrites each term state as a pure- 
state mixture. □ 

Lemma 3. An uncorrelated state p\ <8> P2 has only trivial twins. 

Proof is an immediate consequence of the fact that the tensor factors of a nonzero 
uncorrelated vector, say a £g> b, are unique up to an arbitrary nonzero complex number 
a, but if a is replaced by aa, b must be replaced by (l/a)b. 

If two observables are twins for an uncorrelated state, then 

Mpi ® P2 = Pi ® A 2 p 2 . 

If Aipi = api, then, applying the above remark to supervectors, one has p 2 = {l/a)A 2 p 2 . 
□ 

Proof of Theorem 3 now immediately follows from Lemma 2 and Lemma 3. Namely, 
the two groups of terms stated in the Theorem, make up the two terms in (3). □ 

Corollary 3. Nontrivial twin events of a separable mixture (9) are necessarily strong 
twin events. 

Corollary 4- If [Ai, A 2 ) are nontrivial twin observables for a separable mixture (9), they 
are strong twin observables (cf Definition 3), and the mixture terms can be grouped into 
as many biorthogonal groups of terms as there are distinct characteristic values of A\ 
in TZ(pi) (generalization of (19a,b)). 

It is known that if a statistical operator and a Hermitian operator commute, then 
the corresponding state can be written as a mixture so that each term-state has a definite 
value of the corresponding observable [TU]. But, for the same statistical operator, there 
are also mixtures violating this. 

To take an example, let us think of an unpolarized mixture of spin-one-half states: 
p = (1/2) J (in the two-dimensional spin factor space). This statistical operator 
commutes with s z , nevertheless one can write down the mixture 

p = (1/2) ( \x,+)(x, + \ + \x,-)(x,-\) = (1/2)/, 
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in which the term-states do not have a definite value of the z-component. 

It is interesting that in the case of a separable mixture with a nontrivial twin ob- 
servable it is necessarily its term-states that have the sharp detectable values of the 
corresponding observable. 

7. States with Maximally Disordered Subsystems 

Now we turn to the example that is, for illustrative purposes, investigated in this study, 
i. e., to states (statistical operators) p in C 2 ® C 2 . We say that p is an MDS state (one 
with maximally disordered subsystems or rather subsystem states) if p\ = (l/2)ii and 
p 2 = (1/2)I 2 . R. and M. Horodecki have shown ^T] that for every MDS state there 
exist unitary subsystem operators U\ and U 2 such that 

3 

(Ux <g> U 2 )p(u{ ® U ] 2 ) = (1/4) (/ ® I + Y, U(Ti ® <Ti) = T, (20) 

i=l 

where cij, i = 1, 2, 3, are the well known Pauli matrices o~ x , a y and a z \ and it is seen from 
their place in the expression if they are meant for the first or for the second spin-one-half 
particle. 

Further, they have shown that the operator T is a statistical operator (a quantum 
state) if and only if the vector t from R3 the components of which appear in (20) is not 
outside the tetrahedron determined as the set of all mixtures of the four pure Bell states: 

|V£> = (l/2) 1/2 ( |+) |+) + I-) I-)), |^o 3 ) = (V2) 1/2 ( 1+) I") - I") 1+)), (21) 

where | — ) are the spin-up and the spin-down state vectors respectively. 

It is straightforward to see that the three nonsinglet Bell states | tp s ), s = 1,2,3, 
when written in the form (20), are given by t s = —1, and the other two components 
of t equal to +1. The singlet state | V'o) is i n the form (20) determined by all three 
components of t being equal to —1. 

It is also easy to see that for all mixtures one has 

— 1 < *i < +1 £ = 1,2, 3. 

This is a necessary, but not a sufficient condition for T being a state. In other words, 
the tetrahedron is embedded in a cube, in which there are also nonphysical t. In view 
of the LHS of (20), we call T that belong to the tetrahedron: generating MDS states. 

What we want to find out is: Which of the MDS states have nontrivial twins ? For 
those that do have, we want to find the set of all nontrivial pairs of twins. 

It is sufficient to find the generating MDS states T with nontrivial twins, because 
the validity of 

A\T = A 2 T 
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obviously implies 

(u 1 A 1 ut)(u 1 U 2 TUtU$) = (U 2 A 2 UI)(U 1 U 2 TUM), 

i. e., if the generating MDS states have nontrivial twins, then also the generated MDS 
states do have nontrivial twins, and they are immediately obtained. 

As far as the pure generating MDS states (the Bell states) are concerned, the 
first-particle reduced statistical operator p x is equal to (1/2)/!, all nontrivial Hermitian 
operators A\ commute with it, hence they are twins. To evaluate the corresponding 
twin A 2 , one has to read off the antilinear correlation operator U a from (21) having in 
mind (11), and then utilize (13). For the best known Bell state, the singlet state IV'o)) 
e. g., U a takes | +) into |-), and | -) into (- |+)) (cf (21)). If 

Ai = a++|+)(+| + « — | — ><— | +a+-|+)H +(«+_)* |-)(+|, 

o ++ ,o G R, o+_ G C, 

then the twin A 2 has the form: 

A 2 = ct__ |+)(+| + a++\-)(-\ -a+-\+)(- | -(«+_)* |-)(+| . 

Now we turn to the mixtures of Bell states in our search for nontrivial twins. 

8. Mixtures of Bell states 

Viewing statistical operators as super vectors, and utilizing (redundantly, but for the 
sake of better overview) the ket notation for super state vectors (i. e., Hilbert-Schmidt 
operators as normalized super vectors), one can rewrite the generating vectors T given 
by (20) as a biorthogonal expansion with positive expansion coefficients: 

1 711711-% = (l + I>'r 1/2 ( l(l/2) 1/2 /)i® \{1/2)^ 2 I) 2 + 

i=l 

J2\ti\ KV^/V,)!® |^(t i )(l/2) 1 /V < > 2 ) (22) 

i=l 

("sg" denotes the sign), i. e., as a (super state vector) Hermitian Schmidt decomposition. 

One can read off (22) the following canonical entities of the super state vector 
|T||T||-% (cf (10a), (10b) and (11)): 

The first-subsystem reduced statistical super operator p\ has the characteristic 
super state vectors {| (1/2) 1 / 2 / ) 1; | (l/2) l ' 2 a i ) 1 : i = 1,2,3}; the second 
subsystem reduced statistical super operator p 2 has the characteristic state vectors 
{| (l/2) 1 / 2 /) 2 , | sg{t i ){l/2) l / 2 a i ) 2 : i = 1,2,3}; and the common spectrum of p x 
and p 2 is {R = (1 + X^ =1 t 2 )-\ Ri = Rot 2 : i = 1,2,3}. Finally, the antiunitary 
correlation super operator U a maps the enumerated characteristic state vectors of p\ 
into the correspondingly ordered ones of p 2 . 
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9. Nontrivial MDS twins 

Every super operator A\ that commutes with pi, i. e., for which every characteristic 
subspace of the latter is invariant (and no other super operator ), has a twin super 
operator A 2 [HJ- But we are interested only in those pairs (Ai, A2) in which both super 
operators are, what may be called, multiplicative ones, i. e., which have the form 

Axp 12 = Atpxz, A 2 p 12 = A 2 pi2, 

where A p , p = 1,2, are ordinary (subsystem) operators. It is easy to see that a 
multiplicative super operator is Hermitian (in the Hilbert-Schmidt space of supervectors) 
if so is the ordinary operator (in the usual sense) that determines it. 

The basic result of the expounded illustration is given in the following two theorems: 

Theorem 5. Mixed generating MDS states have nontrivial twins if and only if they are 
mixtures of two Bell states (binary mixtures). 

Theorem 6. A) Let us take a binary mixture of two Bell states both distinct from the 
singlet one, and let T, =| ipi){ipi | (cf (21)) be the nonsinglet Bell state that does not 
participate in the mixture. Then the nontrivial twins are: 

At = ah + (3a\ l \ A 2 = al 2 + (3af\ a,j3eR, (3^0, (23) 

where the suffix on o~i refers to the corresponding tensor factor space. 

B) In case of a binary mixture of the singlet state with another Bell state, say 
Ti ee I ipij (ipi I (cf (21)), the twins are: 

At ee ah + /3<tP, A 2 ee al 2 - (3af\ a,p G R, 0^0. (24) 

Proof of the two theorems and of some subsidiary results is given in the Appendix. 
The proof of Theorem 6 that is first given in the Appendix is only of methodological 
significance: it illustrates a method how to evaluate nontrivial twins. In our case of 
binary mixtures T^ 2 \ another method gives a simpler evaluation. It is given at the end 
of the Appendix. 

It is known that any Bell state can be converted into any other one by local unitary 
transformation |I2j. Hence, for proving the existence of nontrivial twins it would 
have sufficed to take mixtures of one pair of Bell states. Theorem 6 is, nevertheless, 
more elaborate because the explicit form of the twins depends on which Bell states are 
involved. 

It is known that all binary mixtures of Bell states are nonseparable except those 
with equal weights. The latter, as easily seen, are examples of Theorem 4, e. g., as one 
can easily ascertain making use of (21), one has: 

(1/2) ((!+)(+ 1® l+X+l) + (l-X-l® !-)(-!)) = 
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(1/2) ( I^X^I + I^)^ |). (25) 

The nonseparable binary Bell state mixtures are distillable even in the single copy 
case [IS] • Unfortunately, there is no simple relation between the (investigated) existence 
of nontrivial twins and distillability as seen from the fact that also rank four mixtures 
are distillable (if and only if one of the weights is larger than 1/2), and they do not have 
nontrivial twins. 

In conclusion, I would like to point out that the entangled pure state case j^j, [I] 
is a well explored illustration for the fact that nontrivial twins can exist on account 
of entanglement. The nonseparable binary Bell state mixtures are another simple 
illustration for this fact. One should have in mind that, as it was seen in Lemma 
3, uncorrelated bipartite states do not have nontrivial twins. Separable states can have 
nontrivial twins if and only if biorthogonal grouping of the terms is possible (cf Theorem 
4). Unfortunately, for the time being, we do not have a necessary and sufficient condition 
for the existence of nontrivial twins on account of entanglement (except in the pure state 
case), let alone a way of generating all of them for a given composite-system mixed state 
(except in the pure state case). 

At last, but not at least, a relation between the reported twin investigation Q3J, 
jlj, and P besides this article) and the mainstream research on entanglement (take 
the cited Bennett et al. articles, the article of Vedral et al., and the cited Horodecki 
family articles as examples) is still lacking. But I believe that there exists a connection. 
Further research will, hopefully, uncover it. 



10. Appendix 

Since we are going to prove the theorems making use of (22), first we must be able to 
recognize the binary mixtures on the Horodecki tetrahedron. 



Proposition A.l. One has a binary mixture T^ if and only if precisely one of the three 
\ti\ values in (22) equals 1. 



A) If U 



+1, \t 



< 1 (where the three values {1,2,3} of i are meant 
cyclically), then the mixture is of two Bell states both distinct from the singlet state. 
If Tj is the nonsinglet Bell state that does not participate in the mixture, one has 
ti + 2 = —ti + \. Finally, the binary mixture T^ 2 ' in question is 



rp(2) _ 

B) If ti = —1, \ti+2 



1 - t i+1 )/2 T i+1 + 



ti+ 2 )/2 



T 4 



i+2- 



(Al) 



< 1 (in the cyclic sense), then one deals with a mixture 
of two states: the singlet state and another Bell state Tj. One has t i+ i = t i+2 , and the 
binary mixture in question is 



T (2) 



(l + t i+ i)/2 



T + 



;i-ti+i)/2 T . 



(A2) 



Both in the cases (A) and (B), can be any number in the interval — 1 < t 



i+1 



< 



+1; equivalently, one can have any point on the corresponding border of the Horodecki 
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tetrahedron (the vertices excluded). 

For proof a few subsidiary results are required. 

Lemma A.l. If among the four numbers {1, \U\ : % = 1,2,3} appearing in the form 
(22) of the generating MDS state T there is one distinct from the rest, then T has no 
nontrivial twins. 

Proof. As clearly follows from the above stated spectrum of p\, the mentioned "one 
number distinct from the rest" corresponds to a nondegenerate characteristic value. 
Assuming that A 1 is a twin, it is a multiplicative superoperator reducing in each 
characteristic subspace of p\. (This is equivalent to commutation with pi ) 

a) Let us take the case when |tj| < 1, % — 1,2,3. Then the first characteristic 
value of pi is nondegenerate, and the corresponding characteristic super state vector has 
to be invariant (up to a constant): 

A 1 (l/2) 1 / 2 / 1 =o(l/2) 1 / 2 / 1 , 

i. e., Ai — a, and the twin is trivial. 

b) Let \ti\ for some value of % be distinct from the other three numbers. Then the 
corresponding characteristic super state vector o" i (l/2) 1 / 2 must be invariant (up to a 
constant): 

A 1 a l (l/2) 1 / 2 = oa l (l/2) 1 / 2 , 
which, upon multiplication with <7j from the right, implies A\ — a again. □ 

Corollary A.l. If a generating MDS state T has nontrivial twins, then for at least one 
value of i: \ti\ — l. 

Proof 'is obvious from Lemma A.l. □ 

Lemma A. 2. Expressing a generating MDS state T written in the form (22) in terms of 
the statistical weights with respect to the Bell states {T k = \ip k ) (ip k \: k — 0, 1, . . . , 3} (cf 
(6)), one has: 

3 

T = ^w k T k = (1/4) I ®7 + (-wi+w 2 + w 3 -Wo)o- 1 ®o- 1 + (w l -w 2 +w 3 -Wo)o-2®o-2+ 

(wi + w 2 -w 3 - w )a 3 <g> a 3 , (A3) 

where 

3 

VA;: w k e[0,l], k = 0,1, 2, 3; = X - 

fc=0 

Proof is straightforward substituting the Bell states in (22) (cf (21) and beneath it). □ 



Hermitian Schmidt Decomposition and Twins of Mixed States 



20 



Lemma A. 3. If one has = 1, % — 1,2,3, for a generating MDS state T in the form 
(22), then it is a Bell state. 

Proof. Each tj has two sign possibilities; altogether there are 2 3 = 8 possibilities. 
A straightforward analysis of each of these, taking into account Lemma A. 2 and 
J2l=o w k = 1, shows that 4 possibilities do not give states. These are: {sg(ti) — + : % — 

1,2,3}, {H }, { — I — }, and { h}. The remaining four sign possibilities give 

the four Bell states: 

{- + +}: T i; {+-+}: T 2 ; {+ + -}: T 3 ; { }: T . 

□ 

Proof of claim (A) in Proposition A.l. Since it is clear from (A3) that the tj as functions 
of Wk are symmetric (in the sense of the cycle {1,2,3}), it is sufficient to take i — 1. 
Then 

3 

—Wi + w 2 + w 3 — w = 1, and ^ — 1. 

This gives w 2 +w 3 = 1, w\ = w = 0, and t 2 = w 3 — w 2 = —t 3 . Hence, w 2 = (1 — 1 2 )/2 
and w 3 — (1 + 1 2 )/2 as claimed. Since < wi,wo < 1, the claimed intervals for t 2 and 
t 3 follow. □ 

Proof of claim (B) of the Proposition. It runs in full analogy with the proof for case 
(A). □ 

Proof of the main claim of the Proposition. It is easy to see that the proofs of claims 
(A) and (B) of the Proposition go through also for the case when \t i+ i\ or \t i+2 \ equals 
one. Hence, one cannot have \ti\ = 1 for precisely two values of i. If it is so for one 
value, then either it is so for all three values, and one has a pure Bell state, or it is so 
for precisely one value of i, then we have a binary mixture. □ 

Proof of Theorem 6. We now assume that for one value of i, \ti\ = 1, and that the other 
two components of t in (22) are by modulus less than one. Then it is sufficient and 
necessary for an observable A\ that defines a superoperator A\ by multiplication (we 
write this as A\ = (Ai»)) to have a superoperator twin A 2 (that is not necessarily 
multiplicative as Ai) that it reduces in the two-dimensional supervector subspace 
spanned by Ii and af~\ If we write A\ = al\ + Y?j=\ flj^P { a ,(3j GR), and multiply 
with this from the left o~f~\ it turns out that the condition amounts to (3j = 0, j ^ i. 
The symmetrical argument gives the symmetrical result. Thus the multiplicative 
superoperators defined by A\ and, separately, by A 2 do have superoperator twins if 
and only if they are of the form 

A 1 = ah + f3<j? ] , A 2 = 7 / 2 + SaW , (AA) 

where a, /3, 7, 5 e R. 
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The mentioned operators are twins of each other if and only if 

(A 2 .) = UaiA^U- 1 . (A5) 

Now we find out the necessary and sufficient conditions when (A5) is valid for the 
operators given by (A4). Since both sides of (A5) are linear operators, we apply them 
to the basis of supervectors {I 2 , cr? : i — 1, 2, 3}: 

(A 2 .)I 2 = 1 I 2 + 8a? ) - 

(U a {Ai*)Ua l )h = U a (ah + (5a?) = al 2 + sg^a? . 
Thus, we obtain the condition 

7 = at, 8 = sg(ti)/3. 
Utilizing the well known relation 

m=l 

we, further, have 

(A 2 .)af = (7/2 + Saf ] )af = jaf + 8(8^1, + £^ m ^ 2) ); 

m 



(U a (A 1 .)U- 1 )af = sg{t ] )U a {ah+^ ) )a J 
sgitjjUafaa? + (3{8ijh +J2^ijm^)) = 

m 

sg^^iasg^j)^ + f3(8 i:j I 2 -^ie ijm sg(t m )a. 



(2) 



For i = j we obtain the condition j — a, and 8 = sg(ti)/3, and, for j ^ i, in addition: 
8 = —sg(tj)sg(t m )j3. Since i 7^ m 7^ j, we know from the Proposition that, irrespective 
of sg(ti), one has —sg(tj)sg(t m ) = sg(ti). Hence, we actually obtain the condition 
expressed by (23) and (24). □ 
The claim in Theorem 5 that binary mixtures do have nontrivial twins is an 
immediate consequence of Theorem 6. 



The mentioned second, simpler, proof goes as follows: According to Lemma 1, a pair 
of opposite-subsystem observables (A ± ,A 2 ) are twins for a composite-system mixture if 
and only if they are simultaneously twins for each of the pure term states. 

Utilizing (13), it is straightforward to evaluate the twins in the operator basis 

consisting of the four supervectors | — )(— |. But for comparison with the results (23) 
and (24) obtained by the Hermitian Schmidt decomposition method, we do this in a 
little bit more difficult way using the form (22) for the Bell states (see their description 
beneath (22)). 
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We can read off the antiunitary correlation superoperator U a from the mentioned 
form (22) of the Bell state. As it was stated before, every first-subsystem observable 

fti G R, % = 1,2,3), is a twin. The corresponding second- 
subsystem twins for the Bell states are: 

T 2 : A 2 
T 3 : A 2 
To : A 2 

Now, in view of the position of the minus sign in A 2 , evidently, utilizing m ^ i ^ 
j 7^ m m G {1, 2, 3}, and < w < 1, the simultaneous twins are: 

wTj + (1 - w)T m : A 2 = a + foai] 

wT + (1 - w)Ti : A 2 = a - Pi<Ji\ 

and A 2 is, of course, the twin of A\ = a + fiiOi. 
In this way proof of (23) and (24) is obtained. 



= al 2 + (3 l af ) -(3 2 a { 2 2) +(3^; 
= al 2 + (3 1 a[ 2) +(3 2 ai 2) -foa^-, 
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